function y = CL_deltare(x)

 if x >= 0
     c = [0 0 0.884 -0.8964 0.2186 0.0323 3E-05];
     p = [x.^ 6, x.^ 5, x.^ 4, x.^ 3, x.^ 2, x.^ 1, x.^ 0]';
 else
     c = [0 0 -0.884 -0.8964 -0.2186 0.0323 -3E-05];
     p = [x.^ 6, x.^ 5, x.^ 4, x.^ 3, x.^ 2, x.^ 1, x.^ 0]';
 end
 
  y = c * p;
  
end